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ABSTRACT 

Correlated variations in the line and continuum emission from active galactic nuclei (AGN) can be used to 
determine the size and geometry of the broad emission line regions (BLRs). We have spectrophotometrically 
monitored a well-defined sample of 28 Palomar-Green quasars in order to obtain measurements of their BLRs 
and to investigate the relationships between quasar luminosity, central black hole mass, and BLR size in AGN. 
Spectrophotometry was obtained every 1-4 months for 7.5 years, yielding 20-70 observing epochs per object. 
Both the continuum and emission line fluxes of all of the quasars were observed to change during the duration 
of the observing program. Seventeen of the 28 objects were observed with adequate sampling (> 20 independent 
observing epochs) to search for correlated variations between the Balmer emission lines and the continuum flux. 
For each of these 17 objects, a significant correlation was observed, with the Balmer line variations lagging those 
of the continuum by ~100 days (rest frame). Our work increases the available luminosity range for studying the 
size-mass-luminosity relations in AGN by two orders of magnitude and doubles the number of objects suitable 
for such studies. Combining our results with comparable published data available for Seyfert 1 galaxies, we find 
the BLR size scales with the rest- frame 5100 A luminosity as /," 70±{).03 -pj^j^ determination of the scaling of the 
size of the BLR as a function of luminosity is significantly different from those previously published, and suggests 
that the effective ionization parameter in AGN may be a decreasing function of luminosity. We are also able to 
constrain, subject to our assumption that gravity dominates the motions of the BLR gas, the scaling relationship 
between the mass of the central black holes and the luminosity in AGN. We find that the central mass scales with 
5100 A luminosity as M oc /,0 5±o.i -pj^j^ inconsistent with all AGN having optical luminosity that is a constant 
fraction of the Eddington luminosity. 

Subject headings: galaxies: active — quasars: emission lines — quasars: general 



terson et al. 1999), and have presented evidence for Keplerian 
motions of the BLR gas (Peterson & Wandel 1999). 

While Seyfert 1 galaxies have been studied successfully us- 
ing reverberation mapping, few similar studies of the more lu- 
minous AGN - the quasars - have been presented. Those that 
have been published have been limited in their success at de- 
termining the properties of the quasar BLRs. Past attempts 
to monitor quasars spectrophotometrically have generally suf- 
fered from temporal sampling and/or flux calibrations that are 
not sufficient for the determination of the BLR size. Zheng 
and collaborators (Zheng & Burbidge 1986, Zheng et al. 1987, 
and Zheng 1988) reported results of a monitoring program ex- 
ecuted on 30 quasars for several years, with a sampling inter- 
val of about one year They found the emission-line flux to 
change in response to the continuum changes, but because of 
the limited time resolution, only an upper limit of about one 
light-year for the BLR size could be deduced. Several other 
groups (Perez, Penston, & Moles 1989, Korista 1991, Jack- 
son et al. 1992, and Erkens et al. 1995) observed samples 
of quasars for approximately one year sampled with monthly 
observations. Each of these groups reported line variations on 
time scales of a few months, but because of the short duration 
of these programs the data were insufficient for more specific 
conclusions. Wisotzki et al. (1998) monitored the gravitation- 
ally lensed quasar HE 1 104-1805 for five years. Although the 
quasar continuum varied considerably, the emission line fluxes 
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1. INTRODUCTION 

Reverberation mapping, observing the degree and nature of 
the correlation between continuum and emission line flux vari- 
ations, is one of the major tools for studying the distribution 
and kinematics of the gas in the broad line region (BLR) of ac- 
tive galactic nuclei (AGN; see reviews by Peterson 1993 and 
by Netzer & Peterson 1997). About 17 low luminosity AGN 
(Seyfert 1 galaxies) have been successfully monitored and pro- 
duce statistically meaningful results (see Wandel, Peterson, & 
Malkan 1999, and references therein). Best studied among 
these is the Seyfert 1 galaxy, NGC 5548, which was monitored 
from the ground for over eight years, and from space for several 
long periods (Peterson et al. 1999, and references therein). Sev- 
eral other Seyfert Is were observed for periods of order 1 year 
or less, and nine Seyfert Is were studied over a period of eight 
years (Peterson et al. 1998a). The measured time lags between 
the emission lines and the continuum light curves in these ob- 
jects can be interpreted in terms of the delayed response of a 
spatially-extended BLR to a variable, compact source of ion- 
izing radiation. While the observations do not uniquely deter- 
mine the geometry of the BLR, they give its typical size which, 
for Seyfert 1 galaxies, is of the order of light-days to several 
light-weeks (~ lO'^-lO'^ cm). Recent studies have shown that 
the time lags determined in NGC 5548 for different observing 
seasons correlate with the seasonal luminosity of the object (Pe- 
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appeared to remain constant. Researchers have also attempted 
to use spectra obtained with the International Ultraviolet Ex- 
plorer (lUE) and reverberation mapping techniques to study 
AGN. In particular, several studies (Gondhalekar et al. 1986, 
Gondhalekar 1990, Sitko 1990, O'Brien & Gondhalekar 1991, 
and Koratkar et al. 1998) deduced BLR sizes of 0.01-2 light- 
years for AGN. However, these results are controversial due 
to insufficient sampling and/or uncertainty regarding the abso- 
lute flux calibrations of the lUE observations (e.g., Bohlin & 
Grillmair 1988a, 1988b). Even the quasar best studied by lUE, 
3C 273, has yielded disputed results when different researchers 
have analyzed similar lUE data sets. Both O'Brien & Harries 
(1991) and Koratkar & Gaskell (1991a) found a measurable and 
similar lag between continuum and BLR variations, while Ul- 
rich, Courvoisier, & Wamsteker (1993) argue that the line vari- 
ations reported in the earlier studies were only marginally sig- 
nificant. Thus, even in the the best-studied case, past results 
have proven controversial. 

To obtain more definite results on the BLR size in quasars, 
we began monitoring a well-defined sub-sample of 28 quasars 
from the Palomar-Green (PG) sample (Schmidt & Green 1983) 
in 1991 March at the Wise and Steward Observatories. Results 
from the first 1.5 years were presented in Maoz et al. (1994; 
Paper I) where it was shown that most quasars underwent con- 
tinuum variations (at 4800 A rest wavelength) with amplitudes 
of 10% - 70%. Balmer-line variability, correlated with the con- 
tinuum variations, was detected in several objects. The prelim- 
inary data showed that the typical response time of the emis- 
sion lines is < 6 months. Reverberation mapping of such ob- 
jects therefore requires several years, with sampling intervals 
of less than a few months. Further results, based on four years 
of monitoring two quasars from our sample, PG0804H-761 and 
PG0953H-414, were presented in Kaspi et al. (1996a; Paper II). 
The measured time lags between the Balmer lines and the con- 
tinuum variations of the quasars suggested that the BLR size 
grows roughly as the square root of the luminosity of the nu- 
cleus. Wandel et al. (1999) confirmed this result using a larger 
sample of AGN. 

During the period of our spectrophotometric project we also 
monitored several quasar samples photometrically in the B and 
R bands. The broad-band results for radio- and optically- 
selected samples were presented in Netzer et al. (1996) and 
Giveon et al. (1999; hereafter Paper III), respectively. The opti- 
cal sample in Paper III consists of 42 PG quasars, including all 
28 objects discussed in the present paper 

The present paper presents final results from 7.5 years of 
spectroscopic monitoring of our sample. In § 2 we describe 
the sample, observations, data reduction, and present the light 
curves. In § 3 we perform a time series analysis to determine 
the response time of the emission-line flux to continuum vari- 
ations. We estimate the BLR size, continuum luminosity, and 
central mass for each quasar in § 4, and in § 5 we discuss the 
relations among these properties. 

2. SAMPLE, OBSERVATIONS, AND REDUCTION 

The sample, observing technique, and the reduction proce- 
dure, are described in detail in Papers I, II, & III. For the sake 
of completeness, we briefly repeat them here, with emphasis on 
the few procedures which, in the final reduction, were modi- 
fied with respect to previous descriptions. The sample consists 

'*IRAF (Image Reduction and Analysis Facility) is distributed by the National 
cooperative agreement with the National Science Foundation. 



of 28 PG quasars. Objects with a northern declination, B < 16 
mag (based on the magnitudes given in Schmidt & Green 1983), 
redshift z < 0.4, and a bright comparison star within 3'.5 of the 
quasar, were selected. Table [l] presents detailed information 
about the sample. Throughout this paper we will use only the 
first four digits in the PG object name to identify individual ob- 
jects. 

The observations were carried out using the Steward Obser- 
vatory (SO) 2.3m telescope and the Wise Observatory (WO) 
Im telescope. Throughout the years different detectors were 
used. At SO we used the Boiler & Chivens spectrograph with 
a 300 1/mm grating, a 4"5-wide sHt and a TI CCD. In 1992 
September this CCD was replaced with a 1200x800 Loral 
CCD. At WO the Faint Object Specti'oscopic Camera (FOSC; 
Kaspi et al. 1995) was used with a 15"-wide slit and a front- 
illuminated 1024x 1024 TI CCD as a detector In 1994 Febru- 
ary the slit was replaced by a 10"-wide slit and the detector used 
was a 1024x 1024 back-illuminated Tektronix CCD. The typi- 
cal wavelength coverage, at both observatories, was ^ 4000- 
8000 A, with a dispersion of ~3.8 A/pixel, and a spectral reso- 
lution of-- 10 A. 

Spectrophotometric calibration for every quasar was accom- 
plished by rotating the spectrograph slit to the appropriate po- 
sition angle so that a nearby comparison star (columns [11] & 
[12] in Table |l]) was observed simultaneously with the object. 
A wide slit was used to minimize the effects of atmospheric 
dispersion at the non-parallactic position angle. This technique 
provides excellent calibration even during poor weather condi- 
tions, and accuracies of order 1% - 2% can easily be achieved. 

Observations typically consisted of two consecutive expo- 
sures of the quasar/star pair. Total exposure times were usually 
40 minutes at the SO 2.3m telescope and 2 hours at the WO Im 
telescope. The spectroscopic data were reduced using standard 
IRAF^ routines. The consecutive quasar/star flux ratios were 
compared to test for systematic errors in the observations and 
to clean cosmic rays. The ratio usually reproduced to 0.5-2.5% 
at all wavelengths and observations with ratios larger than 5% 
were discarded. We verified that none of our comparison stars 
are variable to within ^^2% by means of differential photometry 
with other stars in each field (Paper III), Spectra were calibrated 
to an absolute flux scale using observations of spectrophoto- 
metric standard stars on one or more epochs. The absolute flux 
calibration has an uncertainty of ^^10%, which is not shown in 
the error bars in our light curves. The error bars reflect only the 
measurement and differential uncertainties, which are of order 
l%-2%. 

For each quasar we used all the spectra to produce an average 
spectrum and a root-mean-square (rms) spectrum, defined as: 

where the sum is taken over the spectra, and /(A) is the av- 
erage spectrum (Peterson et al. 1998a). The mean spectra of 
all the quasars in our sample are shown in Fig. ^ An exam- 
ple of a typical rms spectrum and comparison with the mean 
spectrum is shown in Fig. || for PG0844. We used the aver- 
age and rms spectra to choose line-free spectral bands suitable 
for setting the continuum underlying the emission lines, and 
the wavelength limits for integrating the line fluxes. The spec- 
Optical Astronomy Observatories, which are operated by AURA, Inc., under 
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Table 1 
Sample Properties 



Object 
(1) 


RA (1950.0) Dec 
(2) (3) 


z 

(4) 


mg 
(5) 


Mb 
(6) 


Ab 
(7) 


(8) 


N 
(9) 


A (5 100 A) 
(10) 


(11) 


PA., 
(12) 


PG 0026+ 129 


00 26 38.1 


12 59 30 


0.142 


15.3 


-23.7 


0.13 


38 


56 


26.9 ±4.0 


95 


42.0 


PG 0052+251 


00 52 11.1 


25 09 24 


0.155 


14.7 


-24.1 


0.12 


41 


56 


20.7 ±3.7 


92 


153.4 


PG 0804+761 


08 04 35.4 


76 1 1 32 


0.100 


14.5 


-23.7 


0.11 


38 


70 


54.8±10.0 


182" 


315.6" 


PG 0844+349 


08 44 33.9 


34 56 09 


0.064 


15.1 


-22.1 


0.08 


31 


49 


37.1 ±3.8 


124 


36.8 


PG 0953+414 


09 53 48.3 


41 29 58 


0.239 


15.4 


-24.5 


0.00 


38 


36 


15.6±2.1 


104 


31.7 


PG 1048+342 


10 48 56.1 


34 15 23 


0.167 


17.7 


-21.9 


0.00 


30 


1 


5.51 ±0.28 


186 


94.7 


PG 1100+772 


1 1 00 27.4 


77 15 09 


0.313 


15.9 


-24.7 


0.06 


46 


9 


7.72 ±0.23 


68 


84.1 


PG 1202+281 


12 02 08.9 


28 10 54 


0.165 


16.2 


-23.0 


0.03 


37 


9 


7.2 ±2.0 


73 


219.8 


PG 1211+143 


12 11 44.8 


14 19 53 


0.085 


14.4 


-23.4 


0.13 


52 


38 


56.6 ±9.2 


162 


352.2 


PG 1226+023 


12 26 33.4 


02 19 42 


0.158 


12.8 


-26.3 


0.00 


26 


39 


213 ±26 


171 


171.2 


PG 1229+204 


12 29 33.1 


20 26 03 


0.064 


15.5 


-21.7 


0.00 


32 


33 


21.5±2.3 


75 


291.5 


PG 1307+085 


13 07 16.2 


08 35 47 


0.155 


15.6 


-23.5 


0.02 


35 


23 


17.9±1.8 


138 


186.5 


PG 1309+355 


13 09 58.5 


35 31 15 


0.184 


15.5 


-24.0 


0.00 


30 


9 


23.6 ±1.4 


167 


74.2 


PG 1322+659 


13 22 08.5 


65 57 25 


0.168 


16.0 


-23.3 


0.05 


27 


8 


11. 60 ±0.74 


104 


277.5 


PG 1351+640 


13 51 46.3 


64 00 28 


0.087 


14.7 


-23.2 


0.03 


40 


30 


51.4±5.1 


109 


196.5 


PG 1354+213 


13 54 11.5 


21 18 30 


0.300 


17.3 


-23.2 


0.04 


25 


3 


11.14±0.25 


125 


133.4 


PG 1402+261 


14 02 59.2 


26 09 51 


0.164 


15.6 


-23.6 


0.00 


27 


6 


18.73 ±0.87 


171 


334.8 


PG 1404+226 


14 04 02.5 


22 38 03 


0.098 


16.5 


-21.8 


0.01 


27 


8 


11. 40 ±0.91 


77 


223.0 


PG 141 1+442 


14 11 50.1 


44 14 12 


0.089 


15.0 


-22.9 


0.00 


34 


24 


37.1 ±3.2 


128 


347.0 


PG 1415+451 


14 15 04.3 


45 09 57 


0.114 


16.3 


-22.3 


0.00 


29 


3 


15.2± 1.8 


189 


95.2 


PG 1426+015 


14 26 33.8 


01 30 27 


0.086 


15.7 


-22.2 


0.12 


33 


20 


46.2 ±7.1 


106 


341.4 


PG 1444+407 


14 44 50.2 


40 47 37 


0.267 


17.4 


-23.0 


0.00 


27 


5 


9.98 ±0.39 


158 


155.4 


PG 1512+370 


15 12 46.9 


37 01 56 


0.371 


16.4 


-24.5 


0.05 


35 


10 


4.82 ±0.76 


63 


284.9 


PG 1613+658 


16 13 36.3 


65 50 38 


0.129 


14.9 


-23.8 


0.04 


42 


48 


34.9 ±4.3 


101 


164.2 


PG 1617+175 


16 17 56.2 


17 31 34 


0.114 


15.5 


-22.8 


0.15 


39 


35 


14.4 ±2.5 


147 


253.0 


PG 1700+518 


17 00 13.4 


51 53 37 


0.292 


15.4 


-25.2 


0.02 


41 


39 


22.0 ±1.5 


104 


183.5 


PG 1704+608 


17 04 03.5 


60 48 31 


0.371 


15.6 


-25.5 


0.00 


39 


25 


16.6 ±2.4 


114 


188.2 


PG 2130+099 


21 30 01.2 


09 55 01 


0.061 


14.7 


-22.4 


0.17 


47 


64 


48.4 ±4.5 


86 


68.2 



Notes. — 

(1) Object name, as it appears in Schmidt & Green (1983). 

(2) -(3) Right Ascension and Declination (1950). 

(4) Redshift. 

(5) -(6) Median (over all epochs) apparent and absolute.^ magnitudes from Paper III. 

(7) Galactic B-band extinction, Ab, from NED. See § [4.2| . 

(8) -(9) Number of photometric and spectrophotometnc observing epochs. 

(10) Average flux density and its rms at rest wavelength '^5100 A in units of 10"'^ ergs s"' cm"- A"' (see Table ^ column (7) for exact 
wavelength intervals) . 

(1 1) -(12) The angular separation in arcseconds and position angle in degrees of the comparison star (see § El). 
"Comparison star at SO: Rcomp = 119", PAcomp = 64.3 degrees. 



Table 2 

Integration Limits for Continuum Bands and Emission Lines 



Object 


Continuum 


H7 


Continuum 


Continuum 


H/3 


Continuum 


Continuum 


Hq 


Continuum 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


PG0026 


4790-4830 


4920-5045 


5070- 


-5100 


5408- 


-5468 


5508- 


-5654 


5818-5868 


6998- 


-7138 


7304-7728 


7730-7790 


PG0052 


4864-4914 


4926-5118 


5418- 


-5488 


5418- 


-5488 


5518- 


-5698 


5908-5978 


7028- 


-7108 


7398- 


-7724 


7725- 


-7800 


PG0804 


4680-4705 


4710-4850 


4880- 


-4905 


5224- 


-5264 


5266- 


-5474 


5598-5630 


6746- 


-6814 


6940- 


-7412 


7472- 


-7520 


PG0844 


4503^535 


4536^690 


4710- 


-4750 


5048- 


-5078 


5098- 


-5248 


5414-5440 


6708- 


-6764 


6848- 


-7120 


7248- 


-7348 


PG0953 


5238-5278 


5284-5450 


5472- 


-5506 


5846- 


-5892 


5906- 


-6086 


6288-6320 














PG 1211 


4525^560 


4625^755 


4780- 


-4810 


5120- 


-5160 


5180- 


-5340 


5490-5530 


6670- 


-6720 


6900- 


-7310 


7480- 


-7540 


PG 1226 


4925^945 


4946-5130 


5135- 


-5165 


5494- 


-5522 


5530- 


-5756 


5888-5928 


7034- 


-7085 


7250- 


-7788 


7789- 


-7808 


PG 1229 


4500-4560 


4560-4690 


4694- 


-4758 


5052- 


-5088 


5094- 


-5244 


5412-5456 


6562- 


-6654 


6838- 


-7148 


7244- 


-7320 


PG 1307 


4870-4915 


4938-5088 


5104- 


-5154 


5458- 


-5488 


5508- 


-5698 


5864-5914 


7050- 


-7100 


7400- 


-7740 


7750- 


-7800 


PG 1351 


4585^645 


4646^781 


4815- 


-4850 


5048- 


-5128 


5196- 


-5352 


5545-5605 


6638- 


-6770 


6898- 


-7322 


7362- 


-7432 


PG 1411 


4570-4600 


4650-4790 


4805- 


-4835 


5155- 


-5185 


5201- 


-5382 


5520-5555 


6700- 


-6760 


6870- 


-7360 


7490- 


-7560 


PG 1426 


4565^605 


4605^775 


4730- 


-4796 


5105- 


-5145 


5170- 


-5360 


5515-5555 


6680- 


-6730 


6930- 


-7270 


7490- 


-7550 


PG 1613 










5300- 


-5350 


5360- 


-5566 


5740-5790 


7000- 


-7050 


7013- 


-7742 


7743- 


-7803 


PG 1617 


4650-4700 


4745^950 


4930- 


-4960 


5250- 


-5295 


5310- 


-5550 


5672-5722 


6898- 


-6968 


7118- 


-7498 


7540- 


-7580 


PG 1700 










6030- 


-6080 


6120- 


-6410 


6520-6570 


7760- 


-8015 


8245- 


-8689 


8690- 


-8730 


PG 1704 


5800-5865 


5900-6025 


6040- 


-6100 


6460- 


-6520 


6575- 


-6735 


7020-7100 














PG2130 


4460-4490 


4533^710 


4710- 


-4735 


5048- 


-5084 


5094- 


-5240 


5394-5440 


6498- 


-6598 


6798- 


-7118 


7246- 


-7318 



Note. — Wavelengths in units of A at the observer's frame. The bold-font ranges are the continuum bands shown in Fig. 
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L PG0844 
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Fig. 1.— 



Mean spectra of PG quasars in the observed frame. 
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Fig. 2. — PG0844 mean spectrum (top) and rms spectrum (bottom). 



tral intervals for the Balmer lines and the continuum bands on 
both sides of each line are given in Table |[ The line and con- 
tinuum fluxes were measured algorithmically for all epochs by 
calculating the mean flux in the continuum bands, and summing 
the flux above a straight line in fx that connects the continuum 
bands straddling the emission line. This process measures the 
total emission in each line, i.e., the flux of the broad component 
of the line together with its narrow component. The uncertainty 
in the line flux was estimated by propagating the uncertainty in 
the setting of the continuum levels, determined from the stan- 
dard deviation of the mean in the continuum bands. To this 
we added in quadrature our estimate for the uncertainty in the 
differential spectrophotometry. 

For the purpose of time series analysis, we generally chose a 
region just blue of H/3 for the continuum band. In several cases 
where the blue side of H/3 has an inferior light curve (due to, 
e.g., atmospheric absorption bands and/or superposed emission 
line blends) we chose continuum bands around H7. Table ^ lists 
in bold font the continuum band chosen for each object. Both 
the continuum bands and the emission lines can be affected by 
the contribution of variable Fe II lines. However, correcting for 
this effect is beyond the scope of the present study. Similarly, 
our Balmer line measurements might be contaminated by nar- 
row forbidden emission-lines. Though we measure Hf3 without 
the [Olll] AA4995, 5007 lines (since there is a clear separation 
between them) we cannot do so for the [Nll]AA6548, 6583 
lines blended with the Ha line. However, we note that con- 
tamination from narrow lines should not affect our time series 
analysis since this analysis relies on the variable parts of the 
lines, and the narrow lines have never been observed to vary in 
AGN over relevant time scales (e.g., Kaspi et al. 1996b). 

The AGN continuum dominates the emission-line flux and 
host galaxy starlight contribution to the broad-band photome- 
try. To improve the sampling of the continuum light curves, 
we therefore merged the photometric B band light curves of Pa- 
per III into the spectrophotometric continuum light curves (see 
Paper I). The photometric campaign lasted 7 years and we ex- 
tended the observations of the present 28 objects for six more 
months to better overlap with the spectrophotometric campaign. 
Those last photometric observations were reduced and added to 
the light curves as explained in Paper III. 



The photometric and spectrophotometric light curves of each 
object were intercalibrated by comparing all pairs of spec- 
trophotometric and photometric observations separated in time 
by less than a few days (typically ^5 days, and in some cases 
up to ^10 days). This produced 10 to 15 pairs of points per ob- 
ject. A linear least-squares fit between the spectrophotometric 
continuum fluxes and the photometric fluxes (the magnitudes of 
Paper III translated to fluxes) was used to merge the two light 
curves. 

For each object, the total number of photometric and spec- 
trophotometric observations is given in Table ^ columns (8) & 
(9) (the sum of the two columns gives the total number of points 
in the continuum light curve). Note that for particular emission- 
lines in a few objects, some data points are missing because of 
cosmic-rays contaminations, insufficient wavelength coverage, 
or low S/N. Light curves for the 17 objects that have more than 
20 spectrophotometric data points are presented in Fig. ^ with 
the data listed in Table ^ and Table Q^. Each of the remaining 
11 objects in our sample have less than 10 spectrophotomet- 
ric observations over the 7.5 year period. Such rare sampling 
is insufficient to determine the line-to-continuum lag and these 
objects will not be discussed further in this paper. 

3. TIME SERIES ANALYSIS 

The continuum variability properties of our quasar sample 
are extensively discussed in Paper III. All 17 objects underwent 
optical continuum variability (at 5100 A rest wavelength) of 
25%-150%, measured relative to the minimum flux. In order 
to compare the line variability to the continuum variability in 
each quasar, we define an intrinsic normalized variability mea- 
sure, = 100vo^-<5^//, where / and a are the average and 
the rms of the flux in a given light curve, respectively, and 5 
is the mean uncertainty of all data points in that light curve. 
This measure is less sensitive to noise and outliers than the 
minimum-to-maximum flux ratio. The two cases where - 6^ 
is slightly negative (H7 in PG 1426 and PG 1617) are consistent 
with no variability. Table ^ summarizes for the 5100 A rest 
frame continua band and all of the available emission-line light 
curves. In general we find that all of the Balmer emission lines 
in a specific object have about the same level of variability. In 



'A machine-readable version of these data is available upon request from the authors, or at URL: ittp://wise-obs. tau.ac.il/~shai/PG/ 
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Fig. 3. — Light curves for 17 PG quasars. Circles are spectrophotometric data from WO, squares are spectrophotometric dataJrom SO, triangles are pho- 
tometric data from WO. Continuum flux densities, ix, are determined from the wavelength bands listed in bold-font in Table and are given in units of 
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all objects, the variability of the emission hnes is smaller than 
the continuum variability. In the few cases where those trends 
do not hold, it appears that they may be masked by noise (e.g., 
most cases in which the emission line varies with a larger am- 
plitude than the continuum involve H7, which has lower S/N 
than Ha and H/3). 

The main goal of this project is to determine the time lag 
between the continuum and the Balmer-line flux variations in 
high-luminosity AGN. To quantify this lag, we use two meth- 
ods for correlating the line and continuum light curves. The first 
method is the interpolated cross-correlation function (ICCF) of 
Gaskell & Sparke (1986) and Gaskell & Peterson (1987), as 
implemented by White & Peterson (1994; see also a review by 
Gaskell 1994). The second method is the z-transformed dis- 
crete correlation function (ZDCF) of Alexander (1997) which 
is an improvement on the discrete correlation function (DCF; 
Edelson & Krolik 1988). The ZDCF applies Fisher's z transfor- 
mation to the correlation coefficients, and uses equal population 
bins rather than the equal time bins used in the DCF. The two 
independent methods are in excellent agreement for these data 
and we will use only the ICCF results in the following analyses. 

The uncertainties in the cross correlation lag determination 
can be estimated using the model-independent FR/RSS Monte 
Carlo method described by Peterson et al. (1998b). In this 
method, each Monte Carlo simulation is composed of two parts: 
The first is a "random subset selection" (RSS) procedure which 
consists of randomly drawing, with replacement, from a light 
curve of points a new sample of points. After the points 
are selected, the redundant selections are removed from the 
sample such that the temporal order of the remaining points 
is preserved. This procedure reduces the number of points in 
each light curve by a factor of ^ 1 /e and accounts for the ef- 
fect that individual data points have on the cross-correlation. 
The second part is "flux randomization" (FR) in which the ob- 
served fluxes are altered by random Gaussian deviates scaled 
to the uncertainty ascribed to each point. This procedure sim- 
ulates the effect of measurement uncertainties. Peterson et al. 
(1998b) demonstrate that under a wide variety of realistic con- 
ditions, the combined FR/RSS procedure yields conservative 
uncertainties compared to the true situation. 

All emission-line light curves presented in Fig. | were cross 
correlated with their corresponding continuum ight curves. 
The results of the cross-correlation analysis are presented in 
Table ^. The object name is given in Column (1) and the spe- 
cific Balmer line is listed in column (2). For the ICCF we list 
the centroid time-lag. Teem, (computed from all points within 
80% of the peak value, r„^,x ^) in column (3). We define the 
ICCF peak as the point of maximum correlation and list the 
peak position, Tpeak, in column (4), and the peak value, r,„„x, in 
column (5). In column (6) we list the total number of points in 
the emission-line light curve (the total number of points in the 
continuum light curve is given in Table [l| columns [8] & [9]). 
The uncertainties given for r^,,, and Tpeak were computed with 
the FR/RSS method using ^ 10,000 Monte Carlo realizations 
which were used to build up a cross-correlation peak distribu- 
tion (CCPD; Maoz & Netzer 1989). The range of uncertainties 
contains 68% of the Monte Carlo realizations in the CCPD and 
thus would correspond to Icr uncertainties for a normal distribu- 
tion. In addition to the cross-correlation results, we also list in 
Table || velocity measurements and the average of the observed 
flux for each line. Column (7) lists the mean full- width at half 

^Koratkar & Gaskell (1991a) show the effects of calculating the centroid at 
presented in their study. Therefore, a level chosen at 80% of the peak should be 



maximum (FWHM) of the emission-line, column (8) lists the 
line FWHM measured from the rms spectrum (see § |4.3| be- 
low), and column (9) lists the average of the observed line flux 
and its rms (calculated from the light curves of Table Q). 

Fig. Q presents the cross-correlation functions (CCFs) for the 
40 out of 46 emission-line/continua pairs having peak correla- 
tion coefficients above 0.4 . All 40 CCFs indicate a positive 
time lag of the Balmer lines with respect to the optical contin- 
uum. The time lags are of order a few weeks to a few months, 
and the CCF peaks are highly significant for most lines. In a 
few cases, Tcent and Tpeak are consistent with a zero time lag and 
we attribute this to the fact that these emission-line light curves 
are noisy, and to the conservative error estimate of the FR/RSS 
method. Choosing our best lag determination (see below), i.e., 
avoiding noisy results obtained from weak lines (generally H7), 
we conclude that a time lag has been detected in one or more of 
the Balmer lines for all 17 quasars. 

When comparing the H/3 time-lags for PG 0804 and PG 0953 
obtained in Paper II to the ones deduced here, we find an in- 
crease of about 65% for both objects, although the results are 
consistent at the 2a level. The luminosities of both objects 
also increased during the second half of the monitoring pe- 
riod. Changes in the time lag between different observing sea- 
sons is seen in lower-luminosity AGN (e.g., NGC4151 and 
NGC 5548) and is thought to be a real effect, rather than a re- 
sult of measurement error For instance, the BLR may become 
larger when the object becomes brighter. 



4. SIZE, LUMINOSITY, AND MASS DETERMINATION 



Converting the observable Balmer emission-line time lags, 
line widths, and continuum fluxes into interesting physical pa- 
rameters, namely, BLR sizes, and AGN central masses and lu- 
minosities, is not straightforward. We therefore discuss each in 
some detail below. We analyze our data together with data for 
other AGN. Wandel et al. (1999) have uniformly analyzed the 
reverberation mapping data of 17 Seyfert Is, and deduced time 
lags using the same techniques described in the previous sec- 
tion. Combining their results with the present ones, we obtain, 
for the first time, reliable Size-Mass-Luminosity relations for 
34 AGN spanning over 4 orders of magnitudes in continuum lu- 
minosity. Our analysis is similar to that of Wandel et al. (1999), 
and the few significant differences are noted below. 

We apply below linear regression analysis to the data. Un- 
less otherwise noted, the method is the one described by Press 
et al. (1992), which is based on an iterative process to minimize 
X^, taking into account the uncertainties in both coordinates. In 
some of our data, the uncertainties are asymmetric (e.g., as a 
result of the FR/RSS method and the use of the CCPD), and 
since the fitting method does not account for this, we use the 
mean of the positive and negative uncertainty estimates in a 
given coordinate for each data point. We find that similar results 
are obtained when using the Pearson linear correlation and the 
Spearman rank-order correlation. We present below only Pear- 
son coefficients. 



various levels of the peak value. The ICCFs presented here are smoother than those 
adequate. 
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Table 3 
Continuum Light Curves 

D 

PG 0026 



2448461.5 


22.48 ±0.36 


2448490.6 


20.30 ±0.30 


2448520.3 


19.59 ±0.56 


2448545.8 


20.11 ±0.15 


2448769.6 


21.12± 1.20 


2448836.5 


23.88 ±0.74 


2448869.4 


23.02 ±0.64 


2448890.4 


24.51 ±0.31 


2448925.2 


23.83 ±0.58 


2448949.4 


25.73 ±0.46 



^In units of 10 ergs s ' cm - A ' . 

Note. — Table ^ is presented in its entirety in the elec- 
tronic edition of the Astrophysical Journal. A portion is 
shown here for guidance regarding its form and content. 



Table 4 
Balmer Line Light Curves 



JD 


Flux" 


PG0026-Hc( 




2448836.5 


44.50 ±2.56 


2448890.4 


49.02 ± 1.60 


2448949.4 


60.81 ±1.56 


2448988.6 


53.46 ±0.71 


2449218.6 


51.81 ± 1.57 


2449220.5 


54.04 ±0.99 


2449248.4 


48.81 ±1.03 


2449250.5 


57.43 ± 1.14 


2449278.6 


53.72 ±2.00 


2449281.7 


62.68 ± 0.70 


"In units of 10 '* ergss ' cm ^. 



Note. — Table ^ is presented in its entirety in the elec- 
tronic edition of the Astrophysical Journal. A portion is 
shown here for guidance regarding its form and content. 



Table 5 
Variability Measure o-^f 



Object 


Cont" 


Ha 


H/5 


H7 


PG0026 


14.7 


9.3 


8.0 


8.4 


PG0052 


17.6 


10.1 


11.4 


20.2 


PG0804 


18.2 


5.8 


6.2 


8.7 


PG0844 


10.1 


6.7 


9.7 


10.0 


PG0953 


13.3 




5.5 


6.2 


PG1211 


16.2 


10.3 


12.2 


12.8 


PG 1226 


11.9 


7.7 


7.8 


12.3 


PG 1229 


10.9 


8.1 


12.5 


17.9 


PG 1307 


9.8 


7.9 


12.8 


14.3 


PG1351 


9.8 


4.0 


4.7 




PG1411 


8.5 


3.4 


5.1 


2.8 


PG 1426 


15.3 


7.0 


8.8 




PG1613 


12.4 


5.6 


6.1 




PG1617 


17.3 


7.4 


11.5 




PG 1700 


6.8 




3.2 




PG 1704 


14.5 




10.7 


18.0 


PG2130 


9.3 


6.2 


8.9 


8.9 



"Continuum at rest wavelength ~5I00 A (see Table Q col- 
umn [7] for exact wavelength intervals). 
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Table 6 

Cross-Correlation Results, Velocity Measurements, and Fluxes 



Object 


line 


^cent 


^peak 




Number 


VFWHM(mean) 


VfWHM(rms) 


Flux 






(days) 


(days) 




Obs. 


kms~i 


kms"' 


10~'4ergs s~' cm~2 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


PG 0026 


Ha 


132+^^ 


108!^3 


0.81 


55 


1527 ±99 


1179 ±76 


59.0±5.6 


PG0026 


H/3 






0.84 


53 


2100±140 


1358±91 


12.7±1.1 


PG0052 


Ha 


211!^^ 


1 QO+85 


0.68 


53 


2690 ±5 10 


3280 ±630 


74.6 ±7.9 


PG0052 


H/3 


99i? 


7^+22 


0.75 


56 


3990 ±240 


4550 ±270 


18.0±2.1 


PG0052 


H7 




69!^^ 


0.74 


56 


4050 ±300 


5120 ±380 


8.2± 1.8 


PG 0804 


Ha 


193!?^ 




0.83 


70 


2757 ±55 


2909 ±58 


209 ± 12 


PG0804 


H/3 


151^26 


118^0 


0.85 


70 


2984 ±51 


2430 ±42 


64.5 ±4.2 


PG0804 


H7 




82^^ 


0.74 


67 


3120±150 


2470 ±120 


22.4 ±2.4 


PG 0844 


Ha 


-''^-16 


•"-24 


0.74 


46 


2210 ± 100 


2700 ± 120 


95.3 ±6.5 


PG 0844 


H/3 


1 0+14 a 


1 n+20 
10_i9 


0.79 


48 


2730 ± 120 


2830 ± 120 


30.2±3.0 


PG0844 


H7 


^1 +// 


28!|^ 


0.58 


45 


3210 ±360 


4770 ±530 


13.0±1.9 


PG0953 


H/3 


187+^^ 


207^22 


0.72 


35 


2885 ±65 


2723 ±62 


15.97 ±0.93 


PG0953 


H7 


203!^^ 


213!^^ 


0.54 


34 


3110± 130 


3860 ± 160 


7.12±0.56 


PG1211 


Ha 




AQ+12j 
"^-33 


0.87 


37 


1479 ±74 


1051 ±53 


192 ±20 


PG1211 


H/3 


103!^2 


28!i42 


0.87 


36 


1832±81 


1479 ±66 


54.5 ±6.7 


PG 1211 


H7 


145!^^ 




0.86 


34 


2380 ±150 


2200 ± 140 


27.3 ±3.9 


PG 1226 


Ha 




447+191 


0.71 


34 


2810 ±190 


1330 ±88 


601 ±48 


PG 1226 


H/3 


382ij6^ 


342^Jf 


0.78 


39 


3416±72 


2742 ±58 


172 ±14 


PG 1226 


H7 


3orll 


306!^? 


0.79 


39 


3760 ± 170 


3440 ± 160 


61.1±8.7 


PG 1229 


Ha 


71-46 


40+85 


0.78 


32 


2960 ±110 


3350 ±130 


74.7±6.1 


PG 1229 


H/3 


36!32 


70+1" 


0.80 


33 


3440 ±120 


3490 ±120 


19.8±2.5 


PG1229 


H7 


16!i 


'7+^0 

'-49 


0.81 


29 


3400 ±330 


3000 ±290 


7.8±1.5 


PG1307 


Ha 


179!?45 


00 0+97 
^/«_169 


0.67 


20 


3690 ±140 


3280 ±120 


73.9±6.0 


PG 1307 


up^ 


108!46g 


'3 0+49 

■^^-96 


0.59 


23 


4190±210 


5260 ± 270 


19.2±2.6 


PG1307 


H7 




^-'^-426 


0.64 


23 


4030 ± 240 


3400 ±200 


7.4±1.3 


PG1351 


Ha 


247!jf 


258^^f 


0.49 


29 


1170 ±160 


950 ±130 


76.9 ±3.2 


PG1411 


Ha 


103!^^ 


71+17 

'^-18 


0.87 


24 


2172 ±57 


2135±56 


133.0±4.6 


PG1411 


H/3 


118+^? 


67+1?^ 

"'-31 


0.64 


24 


2456 ± 96 


2740 ± 110 


33 3±1 8 


PG 1426 


Ha 


90!^f 


66!^ 


0.87 


20 


5450 ±150 


4850 ±130 


136.2±9.6 


PG 1426 


H/3 




61+'*i 

"^-66 


0.83 


20 


6250 ±390 


5520 ±340 


30.9 ±2.9 


PG1613 


Ha 


43^0 


T^+185 
■'-'-26 


0.43 


44 


6060 ±210 


1598 ±56 


97.0±5.8 


PG1613 


H/3 


44!20 


29+36 


0.43 


48 


7000 ±380 


2500 ±140 


20.2 ±1.3 


PG 1617 


Ha 


111+31 


83!^5 


0.80 


32 


4270 ±410 


3 160 ±300 


50.4±3.9 


PG1617 


H/3 


78^30 


83^ii 


0.75 


34 


5120 ±850 


3880 ±650 


14.0 ±1.7 


PG1617 


H7 


40+117 


50!if 


0.67 


32 


3 100 ±750 


3110±750 


4.2±1.1 


PG 1700 


H/3 




1 4+296 


0.54 


37 


2180±170 


1970 ±150 


18.88±0.99 


PG 1704 


H/3 


4^7+252 
^•"-391 


^97+263 
-'^'-407 


0.43 


24 


890 ±280 


400 ±120 


3.18±0.41 


PG2130 


Ha 


237+1 


194+|2 


0.59 


64 


2022 ±88 


2236 ±97 


168±11 


PG2130 


H/3 


188+136 


186^144 


0.66 


64 


2410 ±150 


3010±180 


51.2±4.6 


PG2130 


H- 




218+]||i 


0.72 


54 


2850 ± 190 


3050 ±2 10 


22.7±2.7 



''Computed from all points within 60% of the peak value, Ymax, due to the noisiness of the data. 

''The peak of the ICCF measured is the one consistent with the ZDCF results and with the Ha and H7 results. 
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Fig. 4. — Cmss correlation functions: ICCF (solid line) and ZDCF (circles with error bars) for all emission lines with peak correlation coefficients > 0.4. See 
text and Table rj for details. 
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4.1. The Size of the BLR 

Caution must be exercised when using time lags to estimate 
BLR sizes because there are various assumptions and compli- 
cations involved: 

1. The CCF is sensitive to the characteristics of the con- 
tinuum variation, since it is a convolution of the trans- 
fer function and the continuum autocorrelation function 
(Penston 1991; Peterson 1993). Hence, there is a depen- 
dence of the peak and centroid of the CCF on the nature 
of the continuum variability. An illustration of this de- 
pendence is given by Netzer & Maoz (1990). 

2. A real change in BLR size and geometry, over the 7 year 
period, is not out of the question. Using the time-lags 
measured here, and the observed FWHMs, we deduce a 
dynamical (crossing) time of the order of 30 years for the 
more luminous objects and about a year for the lowest 
luminosity Seyfert galaxies. 

3. The measured time lag can be affected by the nonlinear 
response of some emission-line intensities to continuum 
variations. Nonlinear responses are suggested by various 
studies (O'Brien, Goad, & Gondhalekar, 1995; Shields, 
Ferland, & Peterson, 1995) and discussed by Kaspi & 
Netzer (1999). 

4. The ionizing continuum variations may behave differ- 
ently from the observed optical continuum. Due to 
anisotropic continuum emission, the emission line gas 
may be illuminated by an ionizing continuum that is 
poorly represented by the observed optical continuum 
(e.g. Netzer 1987). Furthermore, variable beaming or 
obscuration may be different along the different lines of 
sight from the continuum source to the observer and to 
the gas. Most studies show the continuum varies in phase 
with little or no time delay (albeit, with smaller amplitude 
at longer wavelengths) over a wide wavelength range. 
There are, however, some exceptions (e.g., Peterson et 
al. 1998b) and it is clear that a full understanding of 
the continuum emission processes is yet to be obtained. 
An extreme example is shown by Nandra et al. (1998) 
who found that the X-ray and UV continuum are not 
strongly correlated at zero time lag in NGC 7469, with 
some X-ray emission peaks lagging behind the UV peaks 
by ~4 days. Recently, Maoz, Edelson, & Nandra (2000) 
have demonstrated even weaker correlation at zero lag 
between X-ray and optical variations in NGC 3516. 

5. For a given object, different lines from different ions have 
different time lags. Thus, considering the entire sample, 
the same emission lines should be used in all objects. 
Balmer lines are useful since their lag has been measured 
in most AGN where reverberation mapping has been at- 
tempted. However, as argued by Kaspi & Netzer (1999), 
the theoretical modeling of these lines is very uncertain 
and they are useful mostly because they provide a uni- 
form base for comparing low luminosity Seyferts with 
high luminosity quasars. 

6. There may be flux contributions by other emission lines 
blended with the Balmer lines. These emission lines may 

'when using only H/3 for the analysis, we derive the same results but with s 



have a different time lag and affect the measured emis- 
sion line lag. An example of this potential situation is the 
variable Fe II emission lines near H/3. 

A possible manifestation of some of these difficulties is the 
fact that for some objects, different monitoring campaigns re- 
sult in different time lags. The effect was first noted in the H/3 
time lags found by Peterson et al. (1991; ~7 days) and Netzer 
et al. (1990; -20 days) for the Seyfert galaxy NGC 5548. Other 
examples are the longer time lag deduced for PG0804 and 
PG 0953 in this paper, compared with Paper II (see above), and 
the different time lags in each observing season for NGC 5548 
(Peterson et al. 1999). 

Bearing in mind the above points, we convert the time lags 
(deduced from the ICCF centroid, Tcent from Table ^ of the 
Balmer emission lines directly into BLR sizes (hereafter Rblr) 
after applying a cosmological (1 +z)~' factor We find that time 
lags derived from multiple Balmer lines are consistent with 
each other for individual quasars. The derived linear relations 
are: 

RBLR{Ha) = {\.\9±Q.23)RBiAHI3) + (\3±\9) , (2) 

and 

RBLR{H-i) = {0.96±Q3Q)RBLR{Hf5)-{2>±2>2>) , (3) 

i.e., both are consistent with a slope of 1.0 and a zero intercept. 
In view of this, we adopt the mean time lag measured for Ha 
and H/3 as the best estimate of Rblr '' ■ The relatively low equiv- 
alent width of H7 makes it more difficult to measure and we do 
not include it in the estimate of BLR size. In several objects 
only one of Ha or H/3 is available and we thus use a single line. 
For the Seyfert Is we use the H/3 time lags from Wandel et al. 
(1999) and correct them by the (1 +z)~' factor We have used 
the same method for calculating the time lags as that used by 
Wandel et al. (1999). Values for Rblr for all objects are listed 
in Table ^ column (2). 

4.2. Luminosity Determination 

The exact shape of the spectral energy distribution (SED) of 
most objects in our sample is unknown. Furthermore, much 
of this luminosity is emitted in the unobservable far UV and 
there are still unsolved fundamental issues concerning the shape 
of the continuum (e.g., Zheng et al. 1997, Laor et al. 1997, 
and references therein). Even in those cases where multiwave- 
length data exist, the data are generally obtained at different 
epochs, and variability may render a simple integration over a 
wide wavelength range inappropriate. (See also a discussion 
by Dumont, Collin-Souffrin, & Nazarova; 1998, that demon- 
strates how poorly the SED is determined for even one of the 
best-studied AGN, NGC 5548.) Intrinsic reddening and pos- 
sible orientation effects complicate the situation and it is thus 
unclear what fraction of the total luminosity we are sampling. 

An additional complication is the contribution of the host 
galaxy to the luminosity of the nucleus. While in high- 
luminosity AGN this contribution is negligible, this is not the 
case in low-luminosity objects. Determining the host galaxy 
flux within a given spectrograph aperture is a complicated task, 
attempted only in a few cases (e.g., Peterson et al. 1995). The 
implementation of such corrections is beyond the scope of this 
paper. 

hat larger scatter. 
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Table 7 

Radii, Luminosities, & Masses 



Object 
(1) 



Rblr 
(It-days) 
(2) 



ALa(5100A) 
lO^ergss"^ 
(3) 



M(mean) 
(4) 



M(rms) 
IO^Mq 
(5) 



3C120 
3C 390.3 
Aknl20 

F9 

IC4329A 

Mrk79 

MrkllO 

Mrk335 

Mrk509 

Mrk590 

Mrk817 

NGC3227 

NGC3783 

NGC4051 

NGC4151 

NGC5548 

NGC7469 

PG0026 

PG0052 

PG0804 

PG0844 

PG0953 

PG1211 

PG 1226 

PG 1229 

PG 1307 

PG1351 

PG1411 

PG 1426 

PG1613 

PG1617 

PG 1700 

PG 1704 

PG2130 



40+27 
99 Q+6.3 

37.4!ii 

16.3!?i 
1.4!3.3 

76.7!|^ 
20.0!^-^ 

21 2+2.4 
49+0.6 

113!'? 
134^31 
156!}5 
24.2+17 
I5l!i 

lOl+i 

387!58 

124!^^ 

227!if 

102+38 

95+31 

39!20 
85!^? 

00 + 190 

319!1«4 
200!^^ 



0.73 ±0.13 
0.64 ±0.11 
1.39 ±0.26 
1.37±0.15 
0.164 ±0.021 
0.423 ±0.056 

0.38±0.13 
0.622 ±0.057 

1.47 ±0.25 
0.510±0.096 
0.526 ±0.077 
0.0202 ±0.0011 
0.177±0.015 
0.00525 ±0.00030 
0.0720 ±0.0042 
0.270 ±0.053 
0.553 ±0.016 
7.0±1.0 
6.5±1.1 
6.6±1.2 

I. 72±0.17 

II. 9±1.6 
4.93 ±0.80 

64.4 ±7.7 
0.94±0.10 
5.27 ±0.52 
4.38 ±0.43 
3.25 ±0.28 
4.09 ±0.63 
6.96 ±0.87 
2.37±0.41 

27.1 ±1.9 

35.6 ±5.2 
2.16±0.20 



2 3+'-^ 
^••'-1.1 

34!}i 
18.4^-^ 

8.0!2t 
o.5!li 

5.21-^ 

^•-'"-0.21 

63+"-23 

c 70+0.68 
-'•''*-0.66 

1 7 0+0.44 

4.4!li 

3.9!?i 
94+'' ''2 

^•^^-0.84 

n+013 

1 ca+1.06 
^•-^-'-0.89 

1 2 3+2-3 

- ' •r+0.64 
'-0.65 

5.4!!:? 

22.0!f-^ 

i8.9+;| 

2 16+"-^" 

^•^"-0.83 

18.4+2-8 
4.05!?i? 
55.0!^-^ 

7.5^^ 
28!ji 

4.6!^| 

8.0!3^ 
47+16 
^'-20 
94 1+12.5 

27.3+^-3 
6+13' 

"-13 

3 7+3-1 
14.4!^;^ 



o.65!jj:^^ 0.75; 



3 0+1-9 
■'■^-1.4 

37+12 

■"-14 

18.7!4-« 

8-3!^:i 
0-7!l;| 

io.2!i:^ 
o.77!«:i 

0.38tl^ 

1 Q8+0.34 

n C4+I.O3 

--'•-'^-0.86 

4 Q+2.6 

1 10+1-°^ 

^■^*^-0.98 
14+015 
^■^^-0.09 

i.20!«:^3 
9.4!!:] 

;+0.74 
-0.75 

2 fifi+0.49 
^•""-0.55 

30.2+8-^ 
16.3!1J 

2.7!j:i 

16.4!^i 

^•-'"-0.70 

23.5+3-7 
8-6!^i 
33^ii 
3.0!?;! 

8-8!^i 
37+13 

•"-16 
o 07+1.23 
^•--"-0.88 

15.4!^.v 

*^-'-'-0.81 

20.2!^j 
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Fig. 5. — VfwHMirms) vs. VfwHMii^tean). In this and subsequent figures, PG quasars are denoted by circles and Seyfert Is are denoted by squares. The solid line 
is the best fit to the data. The outlier is PG 1613, which apparently varies only in the core of its very broad lines. The object with the broadest lines is 3C 390.3 . 



Following Wandel et al. (1999), we use XLx at 5100 A (rest 
frame) as our luminosity measure. This luminosity is calculated 
assuming zero cosmological constant, deceleration parameter 
qo = 0.5, and Hubble constant = 75 kms"'Mpc"' (note that 
these values are different from the ones used by Wandel et al. 
1999). This monochromatic luminosity can be translated into 
an integrated 0. 1-1 /im luminosity, assuming a power-law con- 
tinuum of the form cx i/"**-^, as in Papers I & II, by applying 
a correction factor of 3.1. For the 5100 A point, we use the 
continuum measurement on the red side of H/3 (see wavelength 
ranges in column (7) of Table ^). The uncertainty in this value 
is taken to be only due to the variation range of each object, and 
is represented by the rms of the light curve. The observed flux 
has been corrected for Galactic extinction using the Ag values 
from NED^ (see Table |l] column [7]) and a standard extinction 
curve (Savage & Mathis 1979). Monochromatic luminosities, 
ALa(5100 A), are listed in Table |, column (3). 

We have also tried an alternative process to determine the 
optical luminosity. Using line-free bands in the optical spectra, 
we have fitted a continuum to the mean and rms spectra of all 
objects. We integrated those continua over the rest frame wave- 
length of 4500-6000 A to obtain the mean luminosity and its 
range. We have found the results to be very similar to those ob- 
tained by using the monochromatic luminosity. We have also 
examined the consequences of using the H/3 luminosity (rather 
than the optical luminosity) to represent the ionizing luminos- 
ity of each AGN. Here too, we have found the results presented 
below to be mostly unchanged by this choice. Below we use 
only ALa(5100 A) as our luminosity measure. 

4.3. Mass Determination 

To estimate the central masses of the quasars we need to 
assume gravitationally dominated motions of the BLR clouds 

*The NASA/IPAC Extragalactic Database (NED) is operated by the Jet Propul; 
Aeronautics and Space Administration. 



(M w G H'^r). This assumption is presented in Gaskell (1988) 
for the Seyfert 1 galaxy NGC4151 and is also discussed and 
justified by Wandel et al. (1999) and Peterson & Wandel (1999) 
for the Seyfert 1 galaxy NGC 5548 (though see Gaskell 1996 
for a more general approach). There are, however, several com- 
plications regarding the velocity field. Although it is generally 
believed that the line widths represent a Keplerian velocity, it 
is not clear how to convert the observed profile into a velocity 
measure, since the observed profiles are different from the ones 
expected from simple Keplerian orbits of arbitrary inclinations. 
In the absence of such information, it is common to use the rest- 
frame FWHM of the emission line. Because the broad emission 
lines of AGN are composed of a narrow component superposed 
on a broader components, a unique FWHM determination is 
not straightforward. In particular, deblending the components 
is complicated since the line profiles cannot be represented by 
simple analytical functions. Also, the line profiles and widths 
are variable as is evident, for example, in Gilbert et al. (1999) 
and from the spectra obtained in this study. These variations 
can be due to a real change in the BLR velocity field, or other 
reasons, such as blending with weak lines, or observational ef- 
fects (e.g., inaccurate wavelength calibration, or changes in the 
seeing). 

We use two different methods to estimate a typical velocity. 
In the first method, for each quasar we measure the FWHM 
of each line in all the spectra and calculate the mean FWHM, 
VFWHM(inean), and its uncertainty (rms). This empirical way 
of estimating the scatter accounts for all the uncertainties in- 
volved, such as spectral resolution, pixelation of the line pro- 
file, and uncertainties in continuum determination. The values 
are listed in Table ||, column (7). 

The second method, proposed by Peterson et al. (1998a), 
uses the rms spectrum to compute the FWHM of the lines. In 

Laboratory, California Institute of Technology, under contract with the National 




principle, constant features in the mean spectrum (such as nar- 
row forbidden emission Hnes, narrow components of the per- 
mitted emission Hnes, galactic absorptions, and constant con- 
tinuum and broad-line features) are excluded in this method. 
The FWHM measures only the part of the line that varied and 
thus corresponds to the Rblr measured from the reverberation 
mapping'^. For the rms spectrum of each object we have listed 
the emission lines' FWHM, iv//M(rms), in column (8) of Ta- 
ble H There is no simple way to determine the uncertainty of 
VFWHM(rms), and instead we have used the relative uncertainty 
in VfWHMirasan). 

Like the time lags derived using Ha and H/3, the FWHM 
of all lines for a given object are consistent with each other. 
Thus, for the PG quasars we use the average FWHM of Ha 
and H/3. For the Seyfert 1 galaxies we have used the mean and 
rms FWHM quoted in Wandel et al. (1999). The FWHM from 
the rms spectrum versus the mean FWHM from all spectra is 
plotted in Fig. ^. We find the correlation to be highly signifi- 
cant with correlation coefficient of 0.888 . A linear fit gives the 
relation 



5. DISCUSSION 

5.1. Size-Luminosity Relation 

The /?B£R-luminosity relation is presented in Fig. ^. The 
correlation coefficient is 0.827, and its significance level is 
1 .7 X 10"^. A linear fit to the points gives 



Rblr = (32.9!?;^) 



ALa(5100A)\ 
10*4 ergs s-i J 



0.700±0.{)33 



It days (6) 



(solid line plotted in Fig. jg). Considering the Seyfert nuclei 
(log(ALA(5100 A)^ 44.2), or the PG quasars alone, we find 
only marginally significant correlations, probably because of 
the narrow luminosity ranges. A significant correlation emerges 
only when using the whole luminosity range. 

The present result is remarkable for two reasons. First, ear- 
lier studies of this kind found a smaller power-law index in 
the BLR size-AGN luminosity relation (closer to 0.5, e.g., Ko- 
ratkai- & Gaskell 1991b; Wandel et al. 1999). A Hne with this 
slope was fit to the data and is shown as a dashed line in Fig. ^. 
Our result using the combined sample is not consistent with 
previous results. Second, under the assumptions that the shape 
of the ionizing continuum in AGN is independent of L, and 
that all AGN are characterized by the same ionization param- 
eter and BLR density (as suggested by the similar line ratios 
in low- and high-luminosity sources), one expects Rblr oc L^'^ . 
The theoretical prediction is based on the assumption that the 
gas distribution, and hence the mean BLR size, scales with the 
strength of the radiation field. Our present result suggests that 
those assumptions should be re-examined. This is also implied 
from recent models (e.g., Kaspi & Netzer 1999 and references 
therein) which show a wide distribution of properties (such as 
BLR density, column density, and ionization parameter) across 
the BLR of a single active nucleus. Therefore, a range of prop- 
erties may also exist among different AGN, and the above as- 
sumptions of uniform ionization parameter and BLR density for 
all AGN is likely incorrect. If an effective ionization parameter 
can be defined, our result suggests that it may be a decreasing 
function of luminosity. 

'The rms spectra shown by Peterson et al. (1998a) often show narrow forbidden line residuals that are probably due to variable seeing and small shifts in wave- 
length calibration. We see the same effect in our rms spectra as well. It is likely that the same effect operates on the narrow component of the Balmer lines, and 
therefore, the rms spectra are also not completely free of constant narrow components. The degree of this contamination may vary from object to object. 



VFWHM(fms) = (1 .035 ±0.034)vFWHM(mean) - (185 ±93) , 

(4) 

which is plotted as the solid line in Fig. ^ The two velocity es- 
timates are similar and Peterson et al.'s (1998a) preference for 
the FWHM determined from the rms spectra is not empirically 
justified. However, in the subsequent analysis we will use both 
estimates of the broad-line FWHM for each object. 

To determine v we correct VfWHM by a factor of \/3/2, to 
account for velocities in three dimensions and for using half of 
the FWHM. The virial "reverberation" mass is then: 



M= 1.464 X 10^ 



/ ^BLR_\ 

V It days / 



vfwhm 
lO^kms-i 



(5) 



The results of this calculation are listed in columns (4) and (5) 
of Table ^for Vfw'HM(mean) and Vfiy/fM(rms), respectively. 
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XL^(5100A) [ergs s"^] 



Fig. 7. — FWHM — luminosity relations. Top: FWHM derived from the rms spectium of each object. Bottom: FWHM derived from averaging the FWHMs in all 
spectra for a given object. Solid lines are the best fit to the data after excluding the outlier at the lowest luminosity (NGC 4051). 



5.2. FWHM-Luminosity Relation 

Our sample allows us to to re-address the issue of the 
velocity-luminosity relation in AGN. Shuder (1984) noted that 
the FWHM of the Balmer lines increases with luminosity for 
a sample of 25 AGN. Wandel & Yahil (1985) found a correla- 
tion coefficient of 0.5 between the H/3 full width at zero inten- 
sity (FWOI) and the 4000 A continuum luminosity for a litera- 
ture compilation of 94 AGN. Joly et al. (1985; and references 
therein) also reported a weak correlation between the FWOI of 
H/3 and the optical luminosity in a similar collection of objects. 
Boroson & Green (1992) note the FWOI of H/? is quite sensi- 
tive to noise and depends strongly on the quality of the Fe II 
subtraction. These authors found an anticorrelation coefficient 
of -0.275 between the H/3 FWHM and absolute V magnitude, 
significant at the 99% confidence level, in a sample of 87 PG 
quasars. In a recent work Stirpe, Robinson, & Axon (1999) 
measured Ha velocities for 126 AGN and found them to weakly 
correlate with the luminosity. 

In addressing the issue of the velocity-luminosity relation 
we note that our sample is not complete, especially with the in- 
clusion of the heterogeneous Seyfert sample. The results must 
therefore be treated with caution. The FWHMs from the mean 
spectrum are plotted against the luminosity in the bottom panel 
of Fig. ^. There is no significant correlation between the vari- 
ables. However, if we omit the data point for the narrow line 
Seyfert 1, NGC 4051, which clearly deviates from the other 



points and highly influences the fit, we do obtain a marginally 
significant correlation. The correlation coefficient is -0.386, 
and its significance level is 2.7 x 10"^ . A linear fit gives 

v,^HMimean)= (4108!^^) ^^^^i^ ) kms"' . 

\ / \^ lO'*'*ergss"' / 

(7) 

FWHMs from the rms spectrum versus luminosity are plotted 
in the top panel of Fig. y\ The correlation coefficient is -0.540 
has a significance level of 1 .2 x 10"^ and the best linear fit to the 
data is 

, , ,..,..5v,r AL.(5100A) X-'""^"°'°. 

(8) 

Thus, we find significant anticorrelations, in our incomplete 
sample, between the FWHM of the Balmer lines and the lu- 
minosity of the objects such that VfWHM oc /,-0-27±o.oi q^j. 
result indicates a stronger anticorrelation coefficient when 
VFWHM(rms) is used. While the above studies used only one 
epoch for each quasar in their measurements, our study uses the 
means of the luminosities and velocities for each object over 
7.5 years, and averages both Ha and H/3. The PG quasars' 
vpwHM found in this study agrees with the ones found by Boro- 
son & Green's (1992) to within ^10%, except for two objects 
(PG 1351 and PG 1704) for which we measure significantly nar- 
rower lines. If we exclude those two objects, the anticoiTelation 




becomes weaker and its significance decreases. As the correla- 
tion we find is opposite to those found by previous studies, the 
issue of the velocity-luminosity relation needs further investi- 
gation, which is beyond the scope of this paper 

5.3. Mass-Luminosity Relation 

Our two velocity estimates prod uce two mass estimates for 
each object (see Table ^ and § [4.3[ ). Our mass estimates based 
on the determination of the FWHM from the mean spectra are 
plotted versus luminosity in the bottom panel of Fig. ^. The 
correlation coefficient between these two parameters is 0.646 
and has a significance level of 3.7 x 10"^. A linear fit gives 



\ 10^4 ergs s-i J 



0.545±0.036 



Ma 



M(mean)= (5.7 

(9) 

and is plotted as a solid line in the diagram. 

The mass estimates based on the determination of the 
FWHM measured from the rms spectra are plotted in the top 
panel of Fig. ||. We find correlation coefficient of 0.473 with a 
significance level of 4.7 x 10"^. A linear fit to this relation gives 



, ,„7/^ALa(5100A)\ 
M(r^,)= (5.75^.36) X 10 [-j^,,^^) 

and is also plotted as a solid line. 



().402±0.034 



Mq 

(10) 



The results of the above two methods are not consistent. 
Moreover, while it is arguable that u sing the rms spectra to de- 
termine V is a better method (see § 4.3), the mass-luminosity 
correlation based on this measure is less significant. This can 
perhaps be attributed to the fact that the line fluxes in the rms 
spectra are weaker and hence the uncertainty in the correspond- 
ing FWHM larger 

Our M-L relation does not agree with the one found by Ko- 
ratkai- & Gaskell (1991b) of M oc , nor does it agree 

with the one found by Wandel et al. (1999) of M cx L°■^^±o.o^ 
However, when Wandel et al. (1999) use an unweighted linear 
fit they find a slope of 0.54 which is in good agreement with our 
result. We have used a linear regression analysis which takes 
into account the uncertainties in both coordinates (see § 

The fact that the scatter in the mass-luminosity relation is 
larger than that of the size-luminosity and velocity-luminosity 
relations may indicate that luminosity, rather than mass, is the 
variable that mainly determines the BLR size. In fact, from the 
individual size and velocity relations we expect: M oc v^Rblr oc 
^^-o.27-)2^o.7 _ ^0.16^ j ^ weak dependence of mass on lu- 
minosity, with all AGN having similar masses. In practice, we 
have found a somewhat stronger dependence, M (x 5±o i ^ but 
with a large scatter 

The M-L relation can be expressed in terms of the Eddington 
luminosity, LEdd- Roughly estimating the bolometric luminos- 
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ity as Lboi ~ 9ALa(5100A), we obtain an Eddington ratio of 

0.5 

(11) 



Lh„, ^ / ALa(5100A) \ 
LEdd ^ ■ V 1044 ergs s-i J 



The Eddington limit, based on this rough estimate for Li,oi, is 
plotted as a dashed line in Fig. ^. Some of the quasars in our 
sample appear to be approximately at or beyond the Edding- 
ton limit. According to Fig. ^ the Eddington limit traces an 
envelope in the mass-luminosity plane (if we ignore the two 
lower right hand points, PG1700 and PG1704, which have very 
large errors). This depends, somewhat, on the factor chosen 
to estimate the bolometric correction, i.e. the rather uncertain 
slope of the unobserved UV continuum. However, if our cho- 
sen bolometric correction is realistic, we may be seeing a direct 
indication that AGN energy is generated by gas accretion. 

Wandel (1999) reviews three classes of AGN mass estima- 
tion methods. While the X-ray variability method (using the 
shortest time scale for global luminosity variations as the light 
travel time across the Schwarzschild radius) and the accretion 
disk modeling method (deriving an accretion disk model that 
best fits the observed AGN continuum) suggest that the Edding- 
ton ratio increases with luminosity, Wandel (1999) notes that 
the kinematic methods (such as reverberation mapping) have 
yet to show a similar trend. Our reverberation mapping result 
(equation [ll|) indicates for the first time that the Eddington ratio 
increases with luminosity. 

Models that suggest the bulk of the luminosity is due to en- 
ergy release via mechanisms which radiate up to a set fraction 
of the Eddington luminosity are not consistent with the derived 
mass-luminosity relation. One such model is a geometrically 
thin, optically thick, accretion disk that requires L/LEdd < 0.3 
to be self-consistent (Laor & Netzer 1989). Our finding sug- 
gests that the mass accretion rate grows with luminosity much 
faster than the central mass, which would mean very different 
disk properties in low- and high-luminosity sources. Our lu- 
minosity determination is based on the monochromatic flux at 
5100 A. In the thin disk model, this monochromatic flux may 
represents a different fraction of the object's bolometric lumi- 
nosity in AGN of different masses. Thus, it is not clear that the 
thin accretion disk model can be ruled out by the new results. 

6. SUMMARY 

Spectrophotometric monitoring of a large, optically-selected 
quasar sample has shown clear correlations and well-defined 
time lags between the optical continuum and the Balmer- 



line light curves. While the Seyfert 1 galaxies that have 
been studied in this manner all have optical luminosities 
< 1 .5 X lO^"* ergs s~' , the new sample allows us to measure time 
lags in AGN with luminosities up to lO'**' ergs s"'. Our work 
increases the available luminosity range for studying the size- 
mass-luminosity relations in AGN by two orders of magnitude 
and doubles the number of objects suitable for these studies. 
We have combined our results for 17 quasars with data for 17 
Seyfert 1 galaxies having reliable time lag measurements, and 
derived uniform estimates of BLR size, central masses, and lu- 
minosities for the combined sample. 

Our main finding is that the BLR size scales with the 5100 A 
luminosity as ]_P ''0±o.m^ This is significantly different from 
Wandel et al.'s (1999) analysis and is also in contradiction with 
simple theoretical expectations, both suggesting Rblro^ LP'^- 
We have also found that the velocity field of the BLR scales 
inversely with the luminosity, vewhm oc /,-"-27±o.oi Combin- 
ing the measured Rblr with the observed FWHMs, we have 
obtained a mass-luminosity relation for AGN, M oc /,"-5±o i ^ 
which, however, has a large intrinsic scatter. The M-L cor- 
relations are based on two different estimates (mean and rms) 
of the FWHM of the Balmer lines and are not consistent with 
each other, despite the fact that the two measured values for the 
FWHM are generally consistent. Empirically, at least, it is not 
obvious which method of FWHM measurement is preferable. 

Our results show the usefulness of long-term monitoring of 
high-luminosity AGN. There is a need to expand the luminos- 
ity range to include the highest luminosity quasars and this will 
require some 5-10 years of observations. Follow-up studies are 
also needed for some of the results obtained here, such as better 
determinations of the gas distribution in the BLR and the exact 
SED of the quasars. Future work using this sample will include 
the study of time-variable line profiles, lags between continuum 
bands, the intrinsic Baldwin relation, and more. 
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